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We describe a method to construct irreducible baryon operators using all-to-all quark propaga- 
tors. It was demonstrated earlier that a large basis of extended baryon operators on anisotropic, 
quenched lattices can be used to reliably extract the masses of 5 or more excited states in the 
nucleon channel. All-to-all quark propagators are expected to be needed when studying these 
excited states on light, dynamical configurations because contributions from multi-particle states 
are expected to be significant. The dilution method is used to approximate the all-to-all quark 
propagators. Low-lying eigenmodes can also be used if necessary. For efficient computation of 
matrix elements of the interpolating operators, the algorithms should exploit the fact that many 
extended baryon operators can be obtained from the different linear combinations of three-quark 
colour-singlet operators. The sparseness of the diluted noise vectors also afford several computa- 
tion simplifications. Some preliminary results are presented for nucleon effective masses. 
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1. Introduction 

The theoretical determination of the hadron spectrum from first principles is one of the goals 
of the Lattice Hadron Physics Collaboration. Lattice QCD calculations have provided reasonable 
agreement with the physical ground state masses of different baryons (for a review, see Ref. [[T]]) 
and the excited state baryon spectra in the quenched approximation has been reported by several 
groups H, |6[ 0]. In the presence of light, dynamical quarks, it is expected that multi-particle 
states will have to be included for a thorough study of the excited state spectra of hadrons. 

One of the difficulties of studying multi-particle states is the formidable number of contrac- 
tions that one has with conventional point-to-all quark propagators. All-to-all quark propagators 
solve this problem since colour and spin contractions are performed on the sinks and sources sepa- 
rately. One also gains access to disconnected diagrams which are difficult to evaluate with conven- 
tional propagators. In this work, we test the "dilution" method of Ref. [^] for approximating the 
all-to-all quark propagators. 

The low-lying nucleon spectrum obtained through group-theoretical projections ([^]) in the 
quenched approximation has been reported in Refs. [ |Io| , [H] ]. Our aim is to compare the point-to-all 
and all-to-all constructions of these nucleon correlation matrices. We find comparable (or better) 
signals in the single-particle sector and verify that diluted all-to-all propagators can be used to study 
excited baryons. 



2. Construction of Operators/Correlators 

2.1 Point- to- All Construction 

We follow the group-theoretical construction of baryon operators outlined in Ref. . A nu- 
cleon correlation matrix is formed by combining various colour contracted three-quark propagators, 
^(/i"^)(vf^)(T|T) where the indices denote the spin and displacements of the three quarks (sinklsource 
notation), 

r^^)(t\ — 7^(-'') /r'(«"^)(/'P) j_P^{uud){pp) _^{uud)ipp) _^(uud)(pp) 

^ij I'^J - '^MVT'^pvT|'-f(^p)(v|v)(T|T)"^'-'(^|v)(v|7r)(T|T) '^(ai|t)(v| v)(T|7r) v)(v|T)(T|7r) 

_ ^(uud)(pp) _p^{uud)(pp) p^{uud)(pp) ^(uud){pp) \ ,^ 

^(v|v)(T|7r)(M|T) '^(v|7r)(T|v)(M|T) + '^(T|T)(v|v)(Ai|jr) + '^(T|v)(v|T)(M|jr) / ^^■^> 

with the various projection coefficients, vt and cjj^^. Our notation follows closely that of Ref. |^] . 
The (smeared) three-quark propagators are formed from colour contracting the source and sink 
colour indices of the quark propagators, 

G''^^^3Z\a-., ,,-T.(.t)=V^ahce.hr.Q^^^- ,_-(x,f|xo,0)e!? ,^-Jx,t\xo,0) q'^^] ,_y_(jc,f|jco,0) 

{ai\m){Pj\pj]{Yk\Jk)^ ' i-i abc ^aaiplaaip"- ' ' '^bfijp\bP]p^ ' ' ' ' ^cykp\cykp^ ' I ; 

(2.2) 

where Q^^aip\mtp^''^^^'^^^''^^ (smeared) quark propagators as defined in Ref. [||]. 
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2.2 AU-to-AU Construction 

An expression for the nucleon propagator in terms of all-to-all quark propagators is obtained 
by replacing each quark propagator matrix, Q^^\ in the previous expression with 

_ Y N Nm _ 

Qaaip\aa-tpi.xAxo.to) = jim - £^ £ 0[i]ia,>(^'O ® ^[i^ilp ' ^0 ) 

where 7]^^ are the (smeared) diluted random noise sources for source A and (^j^^ are the corre- 
sponding (smeared) solutions (see Ref. [^). The summation over the dilution set is no more than 
A^jjY = NfNcNa in this study and we take N = \ (separately) for each of the three quark propagators. 
(In other words, each of the three quark propagators are estimated independently.) The nucleon 
propagator then takes the form (after some symmetry transformations), 



A B C i j k ^ ^0 



" " " n\ 

<N) / , ,x _ V V V \ \ \ \ V D'i^- 



prW /orO*: I 9#-''' _R-'''*= _r'^^' -R^^-'' -R-''^'" \(r t\ 

'-7m\^^[ABC]pVT"r^'°[CBA]jrvT ^[B/lC]pVT "[ACB]jm ^[CAS]7iVT ^[BCA]jrvT J V-*0) 'Oj 



where 



. N N N N^nKiN^.i _ 
A B C I Ik 



are colour-singlet three-quark operators and 

b\f\x,^)=ijj{x) Uj{x+J)...Uj{x+{p-l)j)d^,,^pj. 

The tilde's indicate that the fields have been smeared. 

The task of computing correlation functions has therefore been reduced to computing the 
source and sink operators, ffj{to;ijk) and ffj{t;ijk) separately and then contracting the remain- 
ing dilution indices to form correlation functions. The latter step does not require large resources 
and hence can be performed as a post-production operation. This scenario is very different from 
the point-to-all construction which requires one to make the correlation functions directly from the 
quark propagators. 

2.3 Optimization 

The dominant computational effort is to make the colour-singlet three-quark operators, fij^gq^y^. 
for each dilution index combination. We are interested in making a correlation matrix with many 
baryon operators. There are three-quark operators that are present in a large number of different 
baryon operators. This is especially the case if we were to make the baryon operators for all differ- 
ent isospins, irreps and rows of irreps at once. Therefore, it is crucial that we only compute each 
three-quark operator once in order to optimize the computational time. In this first study, we shall 
only consider common the three-quark operators that have the same spin indices. Common quark 
displacements will be considered in a subsequent study. 
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3. Simulation/Results 

3.1 Parameters 

We have used the same configurations as in our earUer studies with the conventional, point-to- 
all quark propagators (Ref. [[l^, |lT]]) for comparison purposes. These are the anisotropic = 3) 
quenched, Wilson gauge action (j8 = 6.1) on 12^ x 48 lattices. Anisotropic Wilson fermions with 
pion masses of roughly 700 MeV were used for the valence quarks. 

3.2 Operators 

In this first test, we have focused on the nucleon Gi^ channel with one operator of each "dis- 
placement type" (apart from the single-site case where we have tried all three). These are the 
Single-Site (SS-0,1,2), Singly-Displaced (SD-0), Doubly-Displaced-I (DDI-0), Doubly-Displaced- 
L (DDL-1) and Triply-Displaced-T (TDT-3) type operators. We refer to Ref. for further details 
on these operators. 

3.3 Cost Comparison 

Twenty configurations were analyzed for time-dilution and time-spin-dilution schemes. Time- 
spin-colour-dilution has been tested on a single configuration thus far. Since the traditional method 
requires one to compute the correlation functions directly, we compare the number of quark inver- 
sions required for each computation as an indication of the cost of the simulation time. We note, 
however, that the construction of operators (correlation functions) consumes most of the simulation 
time when using all-to-all (point-to-all) propagators for the heavy quark mass used here. 

Time-dilution The lattice size in the time-direction of the anisotropic lattices used here was 
Nt = 48. The total number of inversions is then 144 = Nt x3 (for three quarks) to generate the 
time-diluted all-to-all quark propagator set. In the conventional point-to-all construction, one needs 
Nc X N(y X 7 = 84 inversions in order to make all of the displaced operators. 

Time-Spin-dilution The number of inversions in this scheme is increased by a factor of 4 to 
576 = 3 X Na X A/f ~ 7 X 84, when we dilute in spin as well as in time. The cost of extra inversions 
for the time-spin-diluted all-to-all propagators is then nearly a factor of 7 compared to the point- 
to-all construction. 

3.4 Effective Masses 

The effective masses are collectively shown in Fig. 1 for time-dilution and Fig. 2 for time- 
spin-dilution. The effective masses with time-dilution alone have large error bars indicating that 
a higher dilution scheme is necessary. On the other hand, the time-spin-dilution scheme shows 
error bars which are already compatible with the point-to-all effective masses on the same config- 
urations. It is important to note that even though the cost of generating extra quark propagators is 
not compensated by a reduction in the error bars, the main purpose of using all-to-all propagators 
is to construct explicit multi-particle operators and correlation-functions on dynamical configura- 
tions, not just to reproduce the quenched spectra of Ref. [ |T0| , pT] |. It is sufficient to have error 
bars which are comparable to those of the point-to-all method for the purpose of analyzing single 
particle states, although further reduction can be achieved by going to a higher dilution scheme 
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t/a, t/a, t/a, t/a, t/a, t/a, 



Figure 1: Time-diluted effective masses for various Figure 2: Time-Spin-diluted effective masses for 
operators in the Gig irrep with isospin 1/2. various operators in the Gig irrep with isospin 1/2. 

if necessary (see Fig. 3 for an example of colour dilution). We also note that the added statistics 
obtained by having source operators on every timeslice on the lattice have been utilized in the 
all-to-all simulations (see Fig. 4). 

3.5 Diagonalization/Excited States 

The main idea behind using a matrix of correlators with many different operators is to find a 
good basis of operators which overlap well with excited states. Since the diagonalization procedure 
involves correlation matrices on two different timeslices, there is a possiblity that the stochastic 
noise on the two timeslices will make the diagonalization unstable. So it is necessary to test whether 
time-diluted noise vectors can be used in the variational method. We have taken one operator of 
each kind (just the first operator for SS) and diagonalized the correlation matrix on timeslice 5. 
The effective masses of the first three states are shown in Fig. 5. Figure 6 shows the corresponding 
results for the point-to-all method. One can see that the diagonalization procedure based on the 
all-to-all method is as stable as the one based on the point-to-all method. The correlation matrix 
for time-dilution alone was too noisy to be diagonalized. This is an example where time-dilution 
alone is not sufficient to even reproduce the point-to-aU results. 

4. Summary 

The efficacy of nucleon operator construction using group-theoretical projections and noise- 
diluted all-to-all quark propagators has been illustrated for selected operators. It was found that 
time and spin dilutions were both necessary to obtain signals for the nucleons that are comparable 
to those of the point-to-all method. We also have a good indication that diluting in colour will 
further reduce the stochastic noise for some of the noisier operators. 

The independent Z4 random sources on each timeslice coming from diluting the stochastic 
source vector did not pose a problem for extracting excited states with time-spin dilution. The 
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Figure 3: The nucleon propagator for time-spin- 
dilution, time-spin-colour-dilution and point-to-all 
quark propagators for the Single-Site-#2 operator in 
the Gig irrep with isospin 1/2 on a single configura- 
tion. 



Figure 4: The error of the nucleon effective mass 
for the SS(0) operator on timeslice 4 as a function of 
the number of sources used on each configuration. 
The dashed line is ^i— . 



signal for the first three states was clearly seen with the limited number of operators and config- 
urations which were used in this first study. The importance of studying the dilution scheme was 
illuminated by the failure of only diluting in time. 

The next step in our study is to construct explicit multi-particle operators with diluted, all- 
to-all propagators. This will enlarge the basis of operators which may be used to obtain a good 
overlap with excited states on dynamical configurations with "light" quark masses. It may also 
become useful to include low-lying eigenmodes together with dilutions (the hybrid method) when 
studying lighter quarks. Work is in progress in these directions. 
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